Langevin simulations of the out-of-equilibrium dynamics of the vortex glass in high-temperature 

superconductors 



Sebastian Bustingorry, 1 Leticia F. Cugliandolo, 2 and Daniel Dommguez 1 

1 Centra Atomico Bariloche, 8400 San Carlos de Bariloche, Rio Negro, Argentina 
2 Universite Pierre et Marie Curie - Paris VI, LPTHE UMR 7589, 4 Place Jussieu, 75252 Paris Cedex 05, France 
and LPT Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex 05 France 

(Dated: February 6, 2008) 

We study the relaxation dynamics of flux lines in dirty high-temperature superconductors using numerical 
simulations of a London-Langevin model of the interacting vortex lines. By analysing the equilibrium dynam- 
ics in the vortex liquid phase we find a dynamic crossover to a glassy non-equilibrium regime. We then focus 
on the out-of-equilibrium dynamics of the vortex glass phase using tools that are common in the study of other 
glassy systems. By monitoring the two-times roughness and dynamic wandering we identify and characterize 
finite-size effects that are similar, though more complex, than the ones found in the stationary roughness of 
clean interface dynamics. The two-times density-density correlation and mean-squared-displacement correla- 
tion age and their temporal scaling follows a multiplicative law similar to the one found at criticality. The linear 
responses also age and the comparison with their associated correlations shows that the equilibrium fluctuation- 
dissipation relation is modified in a simple manner that allows for the identification of an effective temperature 
characterizing the dynamics of the slow modes. The effective temperature is closely related to the vortex liquid- 
vortex glass crossover temperature. Interestingly enough, our study demonstrates that the glassy dynamics in 
the vortex glass is basically identical to the one of a single elastic line in a disordered environment (with the 
same type of scaling though with different parameters). Possible extensions and the experimental relevance of 
these results are also discussed. 

PACS numbers: 74.25.Qt, 64.70.Pf, 61.20.Lc, 74.25.Sv 



I. INTRODUCTION 

Vortex matter in high-temperature superconductors has 
been extensively studied for many yearsi^iii In particular, 
the nature of the dynamics of interacting vortices in their dif- 
ferent phases has attracted much attention due to its relevance 
in technological applications. 

Most studies of vortex systems have focused on the anal- 
ysis of their static properties. It is well known that at low 
temperatures any small amount of disorder destroys the long 
range order of the Abrikosov lattice, 5 leading to different 
glassy phases. For small magnetic field, i.e. small vortex den- 
sity, quasi-long range order is retained in the so-called Bragg 
glass (BG)m2&iL The key characteristics of the BG are the 
logarithmic decay of spatial correlations and the absence of 
dislocations, giving rise to a well defined structure factor. For 
increasing fields, the appearance of dislocation loops signals 
the transformation to a new disordered phase, characterized by 
the rapid destruction of Bragg peaks observed in neutron scat- 
tering. 7 This is the so-called vortex glass (VG)ii2iii Increasing 
the temperature both glassy phases transform into a vortex liq- 
uid (VL). Transport and magnetic measurements showed that 
at low-fields the BG solid melts into the VL through a first 
order transition, while the transition from the VG solid to the 
VL at high fields is continuous. 1213 Whether the latter corre- 
sponds to a true phase transition or a crossover is still under 
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Beyond the nature of the continuous VG-VL transition, im- 
portant efforts were devoted to the interpretation of the irre- 
versibility line (IRL) i 25 ' 26 i 27 Experimentally, the IRL marks 
both the onset of the irreversibility in zero field cooling and 
field cooling magnetization measurements and the separation 



of two transport regimes in the field-temperature diagram. 
Above the IRL, thermal energy dominates, vortices are un- 
able to be pinned, and any amount of current results in a lin- 
ear current-voltage characteristics. Below the IRL, the pin- 
ning energy dominates and flux lines are irreversibly pinned. 
In this case, the resistivity drops to nearly zero with a non- 
linear current-voltage characteristics. Originally, the proposi- 
tion that the IRL signals the continuous VG-VL phase transi- 
tion was supported by experiments showing the scaling of the 
current-voltage characteristic, 28,29 and by simulations in ran- 
domly frustrated 3D XY models without screening showing 
evidence for a finite temperature critical pointil^iSii 3 ^ How- 
ever, contradictions in the experimental scaling of the current- 
voltage characteristics! 2 ^! 3 ^ 2 " on the one hand, and the disap- 
pearance of the finite temperature transition in the randomly 
frustrated 3D XY model when screening is restoredp^i^iii^ 
on the other hand, have raised new questions on how to inter- 
pret the IRL. 

Based on the re-examination of experimental data and the 
results of numerical simulations of the over-damped dynam- 
ics of the London-Langevin model, Reichhart et aM> argued 
that the vortex-glass criticality is arrested at a crossover tem- 
perature. Below this temperature, the relaxation times as ob- 
tained from the study of the resistivity grow very quickly in a 
way that is consistent with the Vogel-Fulcher law commonly 
found in relaxation studies of structural glasses. 33 This study 
suggests that the VL and the VG are then just separated by 
a crossover phenomenon in which relaxation times go be- 
yond the experimental time-window at a glass temperature 
T g (H). Below this crossover line one expects to find an out- 
of-equilibrium system, the VG, with all the properties of more 
standard glasses. 
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In the last decade there has been important progress in the 
understanding of the out-of-equilibrium dynamics of glasses. 
A key characteristic of relaxing glassy systems is the loss 
of stationarity reflected by their aging properties, meaning 
that the system's dynamics depend on the time elapsed af- 
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ter the preparation of the sample 
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As a consequence, 



dynamic correlation functions dependend now on two times, 
the "waiting time", t w , and the time t elapsed during the mea- 
surement. Also the linear response functions of glassy sys- 
tems show aging effects, being dependent on t w and f, and 
they are anomalous in the sense that they are not related 
to their associated correlation functions by the equilibrium 
fluctuation-dissipation theorem (FDT). The relation between 
the two functions, correlation and response, remains, however, 
rather simple. It has been characterized in a number of glassy 
systems, allowing for a kind of classification of the out-of- 
equilibrium dynamics of disordered systems. 35 36 The com- 
parison between linear response and correlation has shed light 
onto the role played by different degrees of freedom in the re- 
laxation of the full system. It allowed for the identification 
of an effective temperature in these non-equilibrium systems 
with slow dynamics. 37 For example, the violation of the FDT 
was studied numerically in driven vortex lattices with random 
pinning^ 

Metastability and aging-like phenomena have been ob- 
served in the nonlinear transport in the VG state of single 
crystal F^S^CaC^Os samples^ffi and analyzed in numer- 
ical simulations iiii^ Also, history dependent effects have 
been found in YBa2Cu3C>7 near the peak effect. 43 Aging 
phenomena has been reported in the magnetization of gran- 
ular Bi2Si"2CaCu208 samples at low fieldsi 4 ^ 5 . However, 
in this case, the effect is due to a chiral glass phase that 
is found in three-dimensional Josephson junctions networks 
with K junctions, 46 thought to model granular high-r c - su- 
perconductors. Theoretically, Nicodemi and Jensen studied 
a simplified two-dimensional lattice Hamiltonian model for 
the effective vortex dynamics to address the issue of equi- 
librium and out-of-equilibrium dynamics. 47 They focused on 
the magnetization relaxation using both one-time and two- 
times correlation functions. They observed that relaxation 
times dramatically increase upon decreasing temperature, 
with a Vogel-Fulcher-like dependence, and established some 
analogies with glass formers and supercooled liquid systems. 
This two-dimensional model emphasizes the effect of random 
quenched disorder but overlooks the importance of the three- 
dimensionality of flux lines and their fluctuations along the 
axis parallel to the magnetic field. Other studies focused on 
the out-of-equilibrium features of vortex creep motion at low 
temperaturesi^Mi 

In the present work, we compare the dynamics of vortex 
matter in high-temperature superconductors with that of struc- 
tural glasses aiming at clarifying the nature of the VG phase. 
We extend and complement our previous study of the out-of- 
equilibrium dynamics of a London-Langevin model. 50 Exten- 
sive numerical simulation of the over-damped dynamics were 
performed to study correlation functions such as the rough- 
ness of the flux lines, their dynamic wandering, the structure 
factor and mean-square displacement all generalized to in- 
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FIG. 1: Sketch of the vortex phase diagram in high-temperature 
superconductors for (a) clean and (b) dirty samples. In (a) the 
Abrikosov lattice phase, prototype of the vortex solid, is located be- 
tween the two critical magnetic fields, H c \ and H c 2, at low tempera- 
tures. At high temperatures the Abrikosov lattice melts into a vortex 
liquid phase at Hm(T). Below H c \ the Meissner phase without vor- 
tex is found. In (b) the Meissner phase at very low magnetic field is 
not shown and the three phases, Bragg glass, vortex glass and vortex 
liquid, are sketched. 



elude two-times dependencies. Focusing on the mean-square 
displacement and its conjugated response function we show 
that the scaling of both functions is multiplicative and we ana- 
lyze how this is related to a non-trivial violation of the FDT. 50 
The outline of the paper is as follows. In Sec. |II| a brief 
description of the vortex phase diagram and glassy proper- 
ties is presented. Section [TO] gives an elaborated description 
of the model and details on the numerical simulations. The 
quantities of interest for the present study are presented in 
Sec. II VI In Sec. [V] we show how the VG phase is approached 
both from the structural and dynamical properties. Finite size 
effects are discussed in Sec. I VII The aging behavior of flux 
lines in the vortex glass and their temperature dependence is 
presented in Sec. l VIII The multiplicative scaling of aging cor- 
relations and responses is described in Sec. lVIIll the violation 
of the FDT, the effective temperature and its relationship with 
the the dynamical arrest of the VL diffusion are discussed. 
Finally, Sec. lIXl is devoted to the discussion of our results. 



II. BACKGROUND 

In this Section we present the general background on which 
we base our study. Our aim is to give here a brief review 
of the vortex matter phases in high-temperature superconduc- 
tors as well as a description of the main features of the non- 
equilibrium relaxation in glassy systems and the violation of 
FDT 



A. Vortex phases in high-temperature superconductors 



ter. 



Several review articles describe the physics of vortex mat- 
Here we briefly summarize the main characteristics 



1.2.3.4 



of vortex phases in high-temperature superconductors, refer- 
ring the reader to these specialized reviews for further 
details. 

In type II superconductors, below a lower critical field 
H c \, the Meissner effect does not allow the magnetic field to 
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penetrate the superconductor sample. Above H c \, the mag- 
netic field penetrates into the superconductor sample, contin- 
uously increasing the flux penetration until the upper critical 
field H C 2 where the sample becomes normal, see Fig.^ The 
magnetic field penetrates the superconductor under the form 
of flux lines (or vortices) each one carrying a flux quantum 
<t>o =hc/2e. Vortices are composed by an inner normal core 
of radius | (the coherent length), surrounded by super-current 
screening over a distance A (the penetration length). In the ab- 
sence of disorder and at low temperatures the vortices form a 
triangular crystal, the so-called Abrikosov latticei 5 ^ 2 ^ due 
to the repulsive vortex-vortex interaction. At high temper- 
ature thermal fluctuations melt the Abrikosov lattice into a 
vortex liquid (VL) phase when the magnetic field Hm < H C 2 
is reached. In the VL vortices move with a non-zero diffu- 
sion constant, leading to a linear current-voltage characteris- 
tics and dissipation. These observations have led to the phase 
diagram in FigQa). The re-entrant behavior of the melting 
line Hm(T) originates in the drastic softening of the elastic 
constants of the Abrikosov lattice when the vortex density is 
lowered close toH c \. 

The high-temperature superconductors reveal new vortex 
physics phenomena related to the relevance of thermal fluc- 
tuations and disorder in these systems. A schematic represen- 
tation of a vortex phase diagram, when disorder effects are 
included, in high-temperature superconductors is presented in 
Fig^b)ii2i£i This type of phase diagram is experimentally ob- 
served in layered systems like YBa2Cu307_ x (YBCO) 5 ^ and 
Bi2Sr2CaCu208+i. 56 In Fig^b) the Meissner phase at very 
low magnetic field is not shown. At high temperature ther- 
mal fluctuations dominate and the VL is observed. At low 
temperature disorder dominates and, when it is short-ranged 
correlated, two glassy phases are observed: the Bragg glass 
(BG)££ at low magnetic field and the vortex glass (VG) at 
high magnetic field. The BG presents a structure factor with 
well defined intensity Bragg peaks at the reciprocal wave vec- 
tors of the triangular lattice. Long-range order, however, is 
destroyed by disorder leading to quasi-long range order with 
spatial correlations decaying logarithmically and power-law 
diverging Bragg peaks*^ At high temperature the BG melts 
into the VL phase through a first-order transition analogously 
to the melting of the Abrikosov lattice. Upon increasing mag- 
netic field or disorder, at low temperatures, the BG transforms 
into the vortex glass (VG) also through a first order transi- 
tion^S 

In the VG phase scenario, disorder destroys long-range or- 
der leading to a short-range amorphous structure. 10,11 While 
the high-temperature VL can be considered to be made of mo- 
bile vortices moving unhindered over the pinning potential, 
the low temperature VG is thought to be an immobile amor- 
phous solid with the flux lines localized on the pinning centers 
in order to minimize the pinning energy. The VG theory pre- 
dicts that collective effects are able to produce infinite energy 
barriers at low temperatures, leading to a strictly zero flux- 
flow resistance with no dissipation. 1011 This leads to interest- 
ing scaling theories for the energies involved in creep motion 
of flux lines £2&Q&L Experimentally, it is very difficult to con- 
clude whether infinite barriers, as predicted by the VG theory, 



really exist (distinguishing between infinite or very high barri- 
ers becomes unfeasible). Moreover, as commonly observed in 
other glassy systems, large energy barriers lead to slow relax- 
ations and below a crossover temperature experiments cannot 
be done in equilibrium. 

B. Aging and violation of FDT 

Aging and the modification of the equilibrium FDT are two 
important properties observed in experiments and simulations 
of glassy systems and also captured in an analytic solutions to 
simple mean-field-like models. 34 - 35 We here briefly describe 
them. 



1. Aging 

Consider the two-times global correlation function 

C(t,t w ) = -^(0(t)0(t w )), (1) 

with 0(t) the global observable of interest, and jV a nor- 
malization factor ensuring C(t,t) — I. In a stationary regime 
the correlation function depends only on the time difference 
Af = t — t w , i.e. C(t,t w ) — C(Af ). Out-of-equilibrium station- 
arity is not ensured. In relaxing glassy systems it is actually 
lost. The correlation then depends on both times, f and t w , 
in such a way that the longer the waiting time the slower the 
decorrelation of the system, i.e. the system is aging. Gener- 
ally, in the long t w limit, the two-times correlation presents a 
separation of time scales, with an initial stationary decay for 
Af <C t w , depending only on Af , followed by the aging regime 
for Af 3> t w . This behavior is usually described by 

C(t,t w )=C st {At)+C ag {t,t w ), (2) 

with C jr (0) = 1 - q EA and C ag (t,t) = q EA , and 

lim Ctf(Ar) = 0, (3) 

lim lim C(t,t w ) = q EA . (4) 

These equations define the Edwards-Anderson order parame- 
ter, q EA . The correlation function first decays from 1 to q EA 
in a time-translation invariant manner for Af <C t w . At longer 
time differences the correlation decays further from q EA to 
depending on both f and t w . This general behavior of the cor- 
relation C is sketched in Fig.|2ja). 

In many cases the aging part of the correlation can be de- 
scribed as C ag (t,t w ) = f [h{t) / h(t w )\, with h{t) a system de- 
pendent monotonic function. This form is found analytically 
in simple coarsening systems below their ordering tempera- 
ture and the p-spin disordered models for fragile glasses. 35 
It also describes very accurately numerical data for Lennard- 
Jones mixtures 62 and light-scattering measurements in col- 
loidal suspensions, 63 among other glass forming particle sys- 
tems. 
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In Fig.|2jb) we present the scaling of the aging part of the 
correlation shown in Fig. [2 a) using a scaling variable t/t w : 



(5) 



The case h{i)<x t" with s any power is usually called "simple 
aging". Note that the short-time stationary regime does not 
scale in Fig. EJb). The type of behavior described by (|5j is 
known as additive aging scaling. 

Sometimes, the simple aging is slightly modified in favor of 
a "sub-aging" or "super-aging" scenario in which h(t) is not 
just a power law of time. In the more common former case 
the ratio h(t)/h(t w ) can be approximated by At/t^ with jx < 1 
in the short time-difference limit. These cases are, of course, 
still compatible with the additive separation in stationary and 
aging regimes described by l|5}- 

A different, less common, functional form describing the 
decay of correlations is given by the multiplicative aging scal- 
ing: 



C(t,t w ) — t w a C a g 



t K 



(6) 



where a > and we have assumed simple aging of the re- 
maining t and t n -dependent function. In this case the scal- 
ing function is multiplied by a f w -dependent correction that 
eventually makes the aging contribution disappear. The scal- 
ing function is such that C ag (t/t w ) ~ f(At)t" for long t w in 
such a way that the correlation reaches, asymptotically, a sta- 
tionary regime. The functional form is found in the auto- 
correlation of the massless scalar field in d > 2j2a systems 
at criticality, 65 ' 66 ' 67 - 68 and a lattice model of a directed poly- 
mer in a 1 + 1 dimensional random environment. 6970 ' 71,72 In 
Figs.|2je) and|2jf) we present a sketch of a correlation func- 
tion with multiplicative aging and its corresponding scaling 
form. 

Note that the difference between the additive and multi- 
plicative scaling is that in the former the correlation function 
in double-logarithmic scale has a well-defined plateau at a 
non-vanishing value separating stationary from aging regimes 
while in the later the correlation eventually becomes com- 
pletely stationary and the aging regime disappears. 

In a diffusive (or anomalous diffusive) problem, another 
type of correlation function is preferred, which is simply the 
displacement 



A(f,f w ) = ([0(f)-0(f w )] 2 ), 



(7) 



that can be rewritten in terms of the correlation C(t,t w ) = 
(0(t)0(t w )) [note that we are omitting here the normalization 
factor jV used in Q] as 



A(t,t w ) = C(t,t)+C(t w ,t w )-2C(t,t w ). 



(8) 



This quantity vanishes at equal times and increases with the 
time difference. If C(t,t) increases with time and is not 
bounded then the growth of A is also unbounded. For this 



type of correlation function it is also possible to write down 
both the additive and multiplicative aging scaling forms: 



A(t,t w ) = A, t (Af)+A a J — 



A(t,t w ) = t%Aa g - 



(9) 
(10) 



in the simple aging case (or with a more general function 
h(t) appearing in A ag in sub or super aging cases). A typi- 
cal displacement has initially a ballistic growth, followed by 
a sub-diffusive regime proportional to At a with a < 1, and 
a normal diffusion regime (a = in the additive aging scal- 
ing scenario). The displacement of a D-dimensional mani- 
fold relaxing in an infinite dimensional embedding space after 
a quench to low-temperature has an additive separation be- 
tween a stationary and an aging regime as in (|5Jl- 73 ' 74 In the 
one-dimensional scalar or TV-component field 64 and Sinai dif- 
fusion 75 - 76 ' 77 the displacement follows instead a multiplicative 
scaling of the form dlOi . In Figs. [2jc)-|2fd) and|2jg)-|2jh) we 
present a sketch of the displacement with additive and multi- 
plicative scaling and their corresponding scaling forms. 



2. The linear response 

In order to measure a linear response the following proto- 
col is commonly used. At t w one applies a perturbation to the 
system by adding the term — hO of intensity h to the total 
Hamiltonian. This perturbation is conjugated to the observ- 
able O. If the perturbation is instantaneous, the linear response 
is defined as 



R(t,t w ) 



_ S(O h (t)) 



8h(t w ) 



(ID 



h=0 



If, instead, the perturbation is held applied during the interval 
[t w ,t] the integrated response associated to the external pertur- 
bation and the observable O is 



X(t,t w ) = 



(O h (t))-(0(t)) 



(12) 



where the sub-index h indicates that O is measured under the 
field. Similarly to what we explained above for the correla- 
tion and displacement, the time-dependence of the integrated 
linear response can be scaled in an additive or a multiplicative 
form: 



(13) 
(14) 



X(t,tw) = Xst{&)+Xag\ - 

. tw 



X{tytw) — t w Xag ( 

Jw 



respectively, where we assumed simple aging. 



3. The fluctuation-dissipation theorem (FDT) 

In equilibrium the integrated response is stationary, 
X{t,tw) = an d m e fluctuation-dissipation theorem 
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FIG. 2: Sketch of the relaxation of correlations and displacements 
and their corresponding scaling forms in glassy regimes, (a)-(b) Cor- 
relation with additive aging scaling; (c)-(d) displacement with addi- 
tive scaling; (e)-(f) correlation with multiplicative aging scaling; (g)- 
(h) displacement with multiplicative scaling. Note that the stationary 
decays at short time-differences in the left panels do not scale in the 
right panels. 



(FDT) implies 



2[1 - C(At )]. Equations 03 and l[T6j then imply 

dA(At) 



dAt 
A (A/) 



2k B TR(At), 
2k B Tx{At). 



(17) 
(18) 



In an aging out-of-equilibrium regime the integrated re- 
sponse also depends on two times, see Eqs. (II 31 and dl 4i . A 
simple generalization of the FDT dl5l reads 



0(r- 



■t w )jr-C{t,t w ) 

at.., 



k B T eS (t,t w )R(t,t w ), (19) 



which gives a definition of a two-times dependent effective 
temperature T e ff. It has been shown that this definition re- 
spects the expected properties of a temperature for systems 
with slow dynamics and a bounded energy. 37 Note that the in- 
tegration of the above relation over a time-interval [t w ,t] can 
take a very complicated form if T e s depends on times in an 
involved manner. If, as it turns out to be in systems with addi- 
tive scaling, it is piecewise constant one finds a linear relation 
between the displacement and the integrated linear response 
in each interval. Systems with multiplicative scaling also al- 
low for the identification of a simple T e ff once the factors t" 
are taken into account. 

The parametric plots %{C) and %(A) have become a useful 
tool to analyze the violation of FDT. 35,36 For systems in equi- 
librium these plots are just straight lines with slope —l/k B T 
and 1 /2k B T, respectively. For glassy systems evolving out 
of equilibrium one has to distinguish those with an additive 
from those with a multiplicative scaling. For the former, the 
parametric plot %(C) [%{A)] shows a -l/k B T [1 /2k B T] slope 
in the quasi-equilibrium short time-difference regime, while 
at long time-differences - for correlation values below the 
plateau at qea - the slope changes to —l/k B T* [l/2k B T*]. 
In the asymptotic, t w — > °° limit the breaking point in x(C) 
occurs at a fixed point [q E A,%EA = \/{k B T){\ - 0ea)]- This 
case is schematically represented in Fig.|5Ja). In the case of 
a system with multiplicative scaling one needs to eliminate 
the factors t~ a in the decaying correlation (t" in the diffus- 
ing displacement) and ?" in the integrated response in order 
to get a stable parametric plot. In other words, tracing %{C) 
[%(A)\ one gets a similar broken line with slopes — l/k B T and 
-\/k B T* [\/2k B T and l/2k B T*] but with a breaking point 
that moves towards (0, l/k B T) [(°°,°°)] with increasing t w . 
This behaviour is schematically represented in Fig.|5Jb). Con- 
structing instead t~ a %(t"C) or, equivalently, t~ a x(t~ a A) a 
stable breaking point is obtained and a parametric plot similar 
to Fig.|3a) is recovered. 



dC{At) 
dAt 
C(Af)-C(0) 



-k B TR(At), 
-k B T X (At), 



(15) 
(16) 



(At > 0) relating correlation and instantaneous or integrated 
responses only through the temperature of the environment. 
In equilibrium the equal-times correlation is constant, and the 
displacement is simply related to the correlation: A(Af) 



III. MODEL 

We model vortices in superconductors as a set of elastic 
lines with an interaction potential Ko(r/X) screened at the 
scale of the London penetration depth A. This model was 
introduced in the theory of flux lattice melting by D. Nel- 
sonr&S it was used to develop the theory of the Bose glass 
phase, 80 and it has been used by several authors in the past 
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FIG. 3: Sketch of the fluctuation-dissipation relation between in- 
tegrated response and correlation: (a) additive scaling with a stable 
breaking point at q£A \ (b) multiplicative scaling with a breaking point 
that drifts towards (0, l/T) in the limit 



and until presentpi2iSiiS2iSi2i There is now a general consen- 
sus in that this model is correct for moderately anisotropic su- 
perconductors like YBCO. Numerical simulations find good 
quantitative agreement with experimental results for this sys- 
tem i 19 i 81 We have chosen to work with such a well-known 
model to build upon previously acquired knowledge and ad- 
dressing questions that have not been considered in the litera- 
ture yet. 



A. Model Hamiltonian 

We consider a model for 3D elastic flux lines in a high- 
temperature superconductor. The model is composed by L 
planes labeled with index z and placed a distance d z apart. The 
direction of the magnetic field B is perpendicular to the planes. 
The magnetic field fixes the vortex density ng = B/<t>o, where 
<t>o = he fie is the flux quantum perpendicular to the planes. 
The system contains N flux lines, where the z'-th flux line 
is characterized by coordinates r,(z), which means that each 
flux line is composed by L elements with two-dimensional in- 
plane coordinates r,(z) = [xi(z),yi{z)]. These elements are 
sometimes called "vortex pancakes". 81 The superconductor 
has anisotropy e = %d%ab = h a b/X c , with and £, c the co- 
herence lengths, and X ab and X c the penetration depths. The 
axes are such that c\\z and ah || r,-. 

Assuming that r,(z) varies slowly with z the Hamiltonian 
for a London model of elastic flux lines is = £ z J^ z , with 

^ = £t/,(Ar fe ) + U d (r iz ) + "£u in {r jz - r iz ), (20) 

i i<j 



where Kq is the modified Bessel function, and with 

£ = {^/AnX ah f. (IT) 
The elastic line energy of the z'-th vortex is 



Ui = -ci{d z rj) 2 dz, 



(23) 



where c/ measures the line tension. A natural choice for d z 
is the distance between CuO planes, in which case one can 
account for the effect of the Josephson coupling 81 using 



U, = ciXj\d z Yi\ for |Ar| > IXj 



(24) 



while keeping the previous expression for |Ar| < IXj, where 
Xj = d z J e is the Josephson length, and Ar,v = r ( - , + i — r,^ mea- 
sures the separation between two adjacent flux line elements. 
The quenched disorder potential due to the impurities is 1 



U d (r) = J d 2 r'u(r')p(\r-r'\ 
where the form factor is 

p{ r )=ie ab /{r 2 +ie ab ), 



and 



{u(r,z)u(r', z '))= Y 8(r-r>)8 Z2 



(25) 



(26) 



(27) 



defines the disorder strength 1 ^ y. 

We model the dynamics of the vortex system with the 
Langevin equation 
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dr iz {t) = 8Jt?[{r u (t)}] , fT 
dt 8r iz 



€(0 . 



(28) 



where rj is the Bardeen-Stephen friction coefficient. The ther- 
mal force ff,(t) satisfies 



(/L(0) = o, 



(29) 



and 



{fizvWhA*")) = 2r lk B T8(t - t'^SffS^,, (30) 
where jU,jLi' = x,y and T is the thermal bath temperature. 



B. Numerical details 



where t//, t/j, and t/,„ are the elastic line energy between two 
pancakes of the same flux line, the interaction of a flux line 
with a quenched disorder potential, and the in-plane interac- 
tion energy between two different flux lines, respectively. Let 
us now describe each term in some detail. 

The repulsive interaction energy between line elements be- 
longing to different flux lines and in the same plane z is ap- 
proximated asiSiiiiH 

U in =1eod z K (r/X ab ), (21) 



There is no standard numerical procedure to compute the 
in-plane interaction (121 1 with periodic boundary conditions. 
Ryu and Stroud worked with a summation calculation over 
image vortices, 81 while Zimanyi and coworkersASiSIi 8 ^ and 
Nordborg and Blatter 82 used instead the periodic extension of 
a simplified Fourier representation. Here we describe the nu- 
merical procedure we used to compute the terms contributing 
to the Hamiltonian HOI . 

The z-planes have dimensions L x x L y and we used peri- 
odic boundary conditions to minimize finite size effects in the 
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planes. To compute the in-plane interaction each direction in 
the z-planes is discretized in N x / y -* points a distance a x u\ apart, 
such that L x /y\ = N x i y ^a x i y ) for the x(y)-direction. The Kq de- 
pendence of the interaction potential, Eqs.(19), corresponds 
to the case of vortices in infinite samples. In a finite sample 
with periodic boundary conditions, the effective potential that 
takes into account the effect of the summation over periodic 
images, Ko(r), has to be calculated. Instead of calculating Kq 
in real space, we found more convenient to work in the recip- 
rocal space, using a discrete Fourier transform of the in-plane 
interaction (12 1 1 that accounts for the periodic boundary con- 
dition. This reads 



U in = 



1 



2eb^3>o 



(31) 



with k and / the indices for the discrete reciprocal space, and 



Ink 



u 



2nl\ 
cos I 

NyJ 



(32) 



The above expressions account explicitly for the periodic 
boundary conditions in L x and L y . Using The numerical value 
of Uin is easily computed in a kl discrete mesh. To calculate 
the real space value of Ui„, a fast Fourier transform FFTW 86 
routine is used. In Fig.|4]we show the Kq and K\ Bessel func- 
tions corresponding to the infinite system, and the effective 
functions Kq and K\ that consider the periodic boundary con- 
ditions, obtained with this method. Different values of X a b are 
shown to follow the general trend. As can be observed, the 
new functions take into account the periodic boundary condi- 
tions: the K\ function, which is proportional to the force be- 
tween vortex pancakes, goes to zero at the edge of the simula- 
tion box (Ljc/2), which is a necessary condition when working 
with periodic boundary conditions. Finally the in-plane force, 
proportional to K\(r/X a b), derived from this potential energy 
is stored in a force table during the computational run. 

Since the z-direction is discretized in L planes, which natu- 
rally correspond to the CuO planes in a model superconductor, 
the elastic interaction is discretized as 



Ui(Ar iz ) = 



for \Ar k \>2Xj, 



for |Ar,v|>2Ay 



(33) 



In order to compute the disorder potential J25I we use the 
same discrete mesh as for the in-plane interaction. The dis- 
crete disorder potential at a given mesh point, indexed by mn, 
is given by 



24^a x a y7 £ £ , 

d= V d L L U mfn'PmV, (34) 

V z m' = \n'=\ 

where u' mn is a random number uniformly distributed between 
and 1 and thus the prefactors enforce Eq. 1271 to hold. The 
discrete form factor is given by 



p m n=Y d Y d [{m-m'fa 2 x + {n-n') 2 a 2 y +2] . (35) 

m'=\ n' = \ 




FIG. 4: Comparison between the analytic Bessel functions Kq (x) and 
Ki (x) and the numerical results using periodic boundary conditions 
and L x /t, a i, = 188. Data is shown up to x = L x /2. 



The indices m and n run over the discrete mesh, i.e. m — 
1,2,... ,N X and n = 1,2,. . . ,N y . In Fig. |5] a typical disorder 
potential realization is shown for two intensities, y = 10~ 4 
and 10~ 5 , later used in the simulations. The force exerted by 
the disorder potential at a given point of the mesh is stored in 
a force table. 

Since the in-plane and disorder potentials are evaluated in 
a discrete mesh but the pancake positions are continuous vari- 
ables in the Langevin dynamics equations, a four-point inter- 
polation method is used. 



C. Parameters of the simulation 

The above model gives a good quantitative descrip- 
tion of the vortex physics of moderately anisotropic high- 
temperature superconductors like Y'&COm^^^ML We there- 
fore choose the values of the parameters to mimic this sys- 
tem: e = 1/5, kab/Zab = 100, Xj/^ ah = 16, and c, = 
£ 2 £o2 [1 + ln(X ab /d z )] /nM The strength of disorder is set to 
y = 10~ 5 , for which case we find that above B cr ~ 0.002H C 2 
the system is in the VG at low temperatures (see SectionlVt. 
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FIG. 5: Quenched disorder potential for two values of the strength 
y. The Abrikosov length on = 26.935 is shown for reference. 

Time is normalized by fo = ^ab 7 ! I £() ' l en g m by the vortex lat- 
tice parameter (or Abrikosov length) ay = [2 < t ) o/(v / 3Z?)] 1 ' /2 , 
energy by £()d z , and temperature by Eod z /kg. We simulate 
N = 56 vortices in a box of size Ioq x 4\/3flo with periodic 
boundary conditions for the in-plane coordinates. The z direc- 
tion is discretized in L = 50 planes with free boundary condi- 
tions. Averages are performed over 10 realizations of disorder. 

IV. QUANTITIES OF INTEREST 

In order to study the out-of-equilibrium dynamics we use a 
two-times protocol. First the system is equilibrated at an ini- 
tial high temperature well inside the VL (7] = 0.3) evolving 
during t = 10 4 steps. Then the system is quenched to a low 
temperature T, where the time count is set to zero. Starting 
with this far from equilibrium initial condition, the system is 
let to evolve during a waiting time t w , after which the quanti- 
ties of interest are measured. 

With the aim of fully characterizing the dynamical proper- 
ties of the vortex system we study several t and f H , dependent 
quantities. First we focus on the study of finite size effects 
by monitoring the roughness of the lines, a quantity that is 
used in the study of interfaces dynamics, 89 here generalized 
to include the t and t w dependence: 

(w 2 (f,f,v)) = ^£([ 5 ^( f )-5x, z (f w )] 2 ), (36) 



where 5x,- z (?) = x, z (f) — jc,-(f) accounts for the displacement 
of the iz line segment relative to the center of mass of the i- 
th elastic line, x,(f ) = L~ l Y lZ Xiz{t)> along one of the two axis 
of the two-dimensional planes on which the pancakes move. 
Another quantity that is useful to analyze finite size effects is 
the dynamic wandering, 90 W(z,?,f„ ), defined as 

W(z,t,t w ) = i£([Ar iZ (f)-Ar fe ( fvv )] 2 ), (37) 



with Ar,,(f ) = r,- z (f ) — r,o(/). This quantity measures how the 
displacement of the vortex pancakes (or line elements) in the 
z-th plane correlates with the bottom plane z = between t and 
t w . Although this was defined for the general vortex model, it 
is worth mentioning that when in-plane interactions and disor- 
der potential are absent, i.e. for the Edwards-Wilkinson (EW) 
equation, 91 the dynamic wandering is related to the height- 
height correlation function studied in surface growth phenom- 



Aging in these systems is best characterized by monitoring 
the mean-square-displacement (MSD) of the pancakes in the 
planes 

B(t,t w ) = ^L(M0 -^-(fw)] 2 ) , (38) 

and the two-times wave-vector dependent density-density cor- 
relation function 

C k (t,t w ) = ± £ (e~* ^ (39) 

that has been analysed in great detail in relaxing glass forming 
liquids. 62 The mean-square displacement B is of the form of 
the generic displacements A discussed and the density-density 
correlation is of the form of the generic decaying correlation 
C discussed in Sec. Ill Bl 

In order to complete the characterization of the aging prop- 
erties and to study the modifications of the FDT in the out-of- 
equilibrium regime, we measure a linear response function by 
applying a random force of the form f IZ = 5s !Z x at time t w on a 
replica of the system, where 8 is the intensity of the perturba- 
tion, and Si z = ±1 with equal probability. 38 - 95 The integrated 
response is 




where xf and x, z correspond to the position evaluated in two 
replicas of the system, with and without the perturbation. The 
equilibrium FDT implies 

B{At)=2k B T X (At), (41) 

where the Af = t — t w argument represents stationarity. In the 
aging regime the FDT is violated and one constructs the mod- 
ified relation 

B(t,t w ) = 2k B T eS (t,t w )x(t,t w ). (42) 

V. VORTEX GLASS 

The vortex phase diagram of a model very similar to 
the one studied in this work was obtained by Zimanyi and 
co-workers, 19 and it has the qualitative behavior shown in 
Fig. |3b). In order to show that we are in the VG region of 
the phase diagram for the present model superconductor, we 
study the behavior of the in-plane vortex structure factor and 
the behavior of dynamic correlations and relaxation times in 
the VL. 




FIG. 6: Structure factor at k\ , the main peak of the Abrikosov lattice, 
and for the disorder intensity y = 1CT 5 as a function of the vortex 
density parametrized by the field B. Insets: a number of superposed 
snapshots of the system. 



A. Structure factor: loss of crystalline order 

In Fig.[6]the vortex structure factor at k\, the principal peak 
of the Abrikosov lattice with parameter ao, is shown while 
varying the magnetic field intensity B (proportional to the vor- 
tex density) and at fixed disorder intensity 7 = 1CP 5 . For this 
case we find that above B cr ~ 0.002i/ C 2 the Bragg peaks disap- 
pear and the flux lines are frozen in a highly amorphous struc- 
ture at low temperatures. Furthermore, B cr is dependent on 
the disorder intensity 7. We therefore choose to study the case 
with B = 0.01// C 2 ^> B cr , which is deep within the VG regime 
at low T for both 7 = 10~ 4 and 10~ 5 . The region with higher 
structure factor S(k\) below B cr presumably corresponds to a 
BG phase, but to be sure one should test the finite size behav- 
ior of the spatial correlations and structure factor. The insets 
clearly show the structural differences between two snapshots 
with different densities. These snapshots are projections of 
the two-dimensional pancake coordinates of the L planes into 
the z = plane. 
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FIG. 7: Correlation Cs(t) for different temperatures and y= 10~ 5 . 
Each curve is labeled by the temperature. The corresponding expo- 
nential fits are shown. 
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FIG. 8: Relaxation times obtained from the exponential fit Cs{t) = 
Ae~'l'' '. The crossover temperature T„, where the VL relaxation time 
departs from the single line behavior t r oc l/T 2 is quoted. The inset 
shows the sudden increase in t r with decreasing temperature. 



B. Diffusion in the VL and dynamical arrest 

Once working at fields B > B cr , upon increasing tempera- 
ture the system should go through the VG-VL crossover line. 
This could be numerically observed as a sudden increase of 
the vortex diffusion. However, since the "pancake" diffusion 
is sub-linear 81 in the VL the data analysis and the identifica- 
tion of the crossover line may be subtle. This can be better 
analyzed starting at high temperatures, from the VL side, and 
studying the relaxation of dynamical correlation functions 

In order to clarify the nature of the dynamic arrest and to 
define a crossover glass temperature T g , we therefore study 
the change in diffusion upon decreasing the temperature in 
the VL phase, by analyzing the MSD defined in Eq. J38I . 



At high temperatures, well inside the VL, the system can 
be easily equilibrated and the MSD depends on Af = t — t w , 
B(t,t w ) = B(At). Two regimes of pancake diffusion can be 
observed in this case. For times longer than a saturation time, 
Af > fjc, diffusion is normal, B(At) ~ At, corresponding to the 
diffusion of the center of mass of each line. At shorter times, 
At < t x , the diffusion of the pancakes corresponds to that of 
line segments in a line in which the longitudinal correlation 
length is still growing. This is the same as considering an 
infinite line L — > °° and it is characterized by sub-diffusion, 81 
B{At)~At l l 2 . 

One aims to extract a characteristic relaxation time from 
the diffusion evolution. To this end, Zimanyi and coworkers 19 
used an exponential fit, C(f) = Ae~ 1 l ty , to describe the corre- 
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lation 



C(f)=e X p{-([r(f)-(r(f))]2)}, (43) 

where r(?) measures the two-dimensional in-plane displace- 
ments of the flux line elements. However, this is only consis- 
tent with a normal diffusion regime, and should be carefully 
interpreted. In their simulations they used an equilibration 
time of 10 s and measured C(t) for the following 10 s time in- 
terval. This is the same as using a waiting time of 10 5 from the 
given initial condition. The exponential fit is valid if t x < 10 5 , 
corresponding to center of mass diffusion. 

As mentioned above, Ryu and Stroud 81 studied the melting 
of the VL using a different model and observed that the pan- 
cake diffusion is sub-linear, B(At) ~ At a with the exponent 
a = 1/2 at high temperatures and a < 1/2 at low tempera- 
tures. This is consistent with our results. 

Within the same scheme, we calculate the evolution of the 
correlation for the largest waiting time used in our simula- 
tions, t w — 10 5 , 

C B (f)=exp{-[B(f,f lv =10 5 )] 2 }, (44) 

and fit it to the form Cfl(f) = Ae^'l'' ', that follows after using 
the sub-linear exponent a = 1/2. In Fig.Qwe show this corre- 
lation in a log-linear representation for different temperatures 
and 7 = 10~ 5 . The fits for Af > 0.6 x 10 5 are also shown. 
From this data we obtain the relaxation time t r presented in 
Fig. [8] The inset shows that there is a sudden increase in 
the relaxation time upon decreasing temperature. For non- 
interacting flux lines the temperature dependence of the relax- 
ation time should be t r oc 1 /T 2 r&^ and the same is expected 
deep in the VL. The main panel in Fig. [8] shows that indeed 
for high temperatures t r grows linearly with 1/T 2 . Below 
a temperature T g w 0.18 we observe that the relaxation time 
t r departs strongly from the 1 /T 2 -dependence and presents a 
remarkable growth. Furthermore, for T < T g is not possible 
to equilibrate the system and correlation functions depend on 
the two times, t and t w , as we will analyze in Sec. lVIII There- 
fore, the temperature T g signals the crossover from the equilib- 
rium high-temperature regime to the out-of-equilibrium low- 
temperature regime. 



VI. FINITE SIZE EFFECTS 

Before analyzing the out-of-equilibrium behavior of differ- 
ent correlation and response functions of the VG, it is neces- 
sary to previously identify finite size effects in the dynamics 
of the system. We here compare finite size effects in the sin- 
gle elastic line dynamics with those in the much more com- 
plex systems we are dealing with, a set of interacting lines 
in the presence of quenched disorder. To this end we study 
both the two-times roughness d36i and the dynamic wandering 
1371 . We also show how finite size effects appear in two other 
quantities of interest, the pancake MSD I38i and the density- 
density correlation function ( 1391 . 
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FIG. 9: (a) Roughness and (b) dynamic wandering of non- 
interacting elastic lines without quenched disorder (EW elastic lines) 
at T = 0.06. Three dynamic regimes are highlighted, corresponding 
to single pancake, single line, and center of mass dynamics. Note that 
different curves with t w = 10 4 , 10 5 and 10 6 (open squares, open dia- 
monds and open triangles, respectively) overlap, showing that these 
quantities are in a stationary regime. 



A. Single lines: the Edwards- Wilkinson equation 

When the in-plane interactions and disorder potential are 
absent from the Hamiltonian J20t . Eq. J28t i reduces to N 
decoupled Edwards-Wilkinson (EW) equations. 91 The EW 
equation was studied in great detail in the context of growing 
interface dynamics. 89 After a transient the roughness of the 
EW interface becomes stationary and it follows the Family- 
Vicsek scaling^ 

<^(Ar)> (45) 

where t x ~ L e is the saturation time, L is the system size, 
and f(x) is a scaling function with f(x) ~ x 2 P for x -C 1 and 
f{x) ~ cons, forx 3> 1. The growth and roughness exponents, 
j3 and £ respectively, are related to the dynamic exponent 9 
through the scaling relation 9 = £//3 . The exponents take val- 
ues £ = 1/2, P = 1/4, and 9=2. This scaling form clearly 
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FIG. 10: Dynamic wandering as a function of the plane index z. 
T = 0.06, L = 50 each curve is labeled by At and the waiting time is 
t w — 10 6 for all. Since W(z,At) saturates at a value proportional to 
the longitudinal correlation length the growth of Si is observed 
for increasing At at fixed z. 



states how the system size is involved in the line dynamics. 
Other quantities, such as the global correlation function, can 
however be non-stationary and show aging. 64 

In order to check how the size and time dependence of our 
elastic lines compare to the EW ones, we set disorder and 
in-plane interactions to zero and we analyse the full At de- 
pendence of the roughness and dynamic wandering once the 
stationary regime has been reached. In Fig. |9j a ) we show the 
time-evolution of the roughness. Three regimes are clear in 
the figure. For Af < f i the roughness grows as (w 2 ) ~ At which 
corresponds to a regime without elastic interactions between 
the line segments. For At > ti the roughness saturates to an In- 
dependent value, while for t\<At< ti the roughness grows as 
(w 2 ) ~ At 1 / 2 . The last two regimes are the ones in the scaling 
function of the Family-Vicsek scaling (145b . where t x = ti. 

To better understand the meaning of a characteristic 
distance-dependent time f- let us analyze the behavior of the 
dynamic wandering in Fig.|9jb). For At <t\ the dynamic wan- 
dering is independent of z and grows as ~ At. This is the 
regime in which the line segment in one plane evolves inde- 
pendently of the line segment in any other plane, i.e. there is 
no truly elastic interaction yet. For At > t\ the line segments 
belonging to two adjacent planes are elastically interacting 
and their relative position saturates to a given value. However, 
at this time, line segments separated by more than a single 
plane are still uncorrelated. This corresponds to the growth 
of a longitudinal correlation length, h\, along the z-direction. 

This correlation grows as §y (At ) ~ Af l l e for At <C t x and satu- 
rates at the system size, 89 §y (Af ) ~ L, for At t x . The time f, 
signals the moment when the correlation length is of the order 
of z, i.e. &i(f z ) ~ z. In Fig.|9jb) the behavior of W(z,t ,t w ) is 
presented for z — 1,3, 10 and z = L = 50, where it is shown 
that the dynamic wandering is essentially a measure of the 
longitudinal correlation length (the short time Af 1 / 2 regime 




FIG. 11: (a) Pancake MSD B and (b) correlation Q for the same 
parameters as in Figs.l9land ll0l ko is the wave vector at the first peak 
of the structure factor. The three dynamic regimes are clear in B, but 
the correlation Ck does not present any signature of them. Symbols 
are as in Fig. [9] 



is not completely developed since the system size is not big 
enough). 

In Fig. ^3 we show the z-dependence of the dynamic wan- 
dering. Since there is no f„,-dependence, we show here differ- 
ent curves labeled by the corresponding Af value for a single 
waiting time t w = 10 6 . For Af < t\, the dynamic wandering 
does not depend on z since line segments are independent of 
each other. For t\ < At < t x two regimes are observed. For 
z < §|| (At) one has W ~ z, while for z > £,\\(At) the dynamic 
wandering saturates at a Af -dependent value. The final growth 
near z = L is due to the free boundary conditions of the line. 
Finally, for Af > t x all the line is correlated and W ~ z for all 
z. 

The growth of the longitudinal correlation length is also 
clear in other quantities. For instance, in Fig. llir a) we show 
the pancake MSD B(t,t w ), where the three dynamic regimes 
are also observed. For Af < t\ independent pancake diffusion 
is observed and diffusion is normal, B(At) ~ Af . For f i < Af < 
t x diffusion is characterized by a sub-linear law B(At) ~ Af 1 / 2 , 
corresponding to the single line regime where every pancake 
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FIG. 12: (a) Roughness and (b) dynamic wandering in the VL at 
T = 0.2 with disorder intensity y = 10~ 4 (cfr. Fig.|5|where the same 
measurements for the EW line are shown). 



elastically feels the others. ' Finally, for At > t x , when the 
correlation length is of the order of the system size, a trend to 
normal diffusion corresponding to the center of mass diffusion 
is observed. 

Lastly, in Fig. II lf b) the density-density correlation function 
Ck (t,t w ) is shown, where Icq is the position of the first maxi- 
mum of the static structure factor. This correlation shows how 
the position of the pancakes decorrelate in time. Note that 
around At ~ 10 4 the pancakes have lost all information about 
the position at t w (At — 0) and Q has effectively decayed to 
zero. There is no signature of the three dynamic regimes in 
the decay of the correlation function. 

In this subsection we demonstrated how the system size is 
involved in different dynamical quantities. We also identified 
three main dynamic regimes. In the following subsections we 
study how these features are reflected in the VL and VG. 



B. The vortex liquid 

Here we return to the case with in-plane interactions and 
disorder, and present results at high temperature, T — 0.2, and 



FIG. 13: (a) Roughness and (b) dynamic wandering in the VG at 
T = 0.05 with disorder intensity y = 10~ 4 (cfr. Figs.l9land ll2t where 
the same measurements in the EW line and the VL, respectively, are 
shown). In (b) the dynamic wandering is shown only for t w = 10 4 . 

with disorder intensity y = 10~ 4 , well inside the VL. The anal- 
ysis of finite size effects in the roughness and dynamic wan- 
dering (Fig. I12> shows that these are similar to the ones in 
the EW line. The values of t z are of the same order. At this 
high temperature and for this disorder intensity one can ar- 
gue that thermal fluctuations are greater than disorder induced 
fluctuations, and the disorder potential in J20i can be basi- 
cally neglected. Moreover, although in-plane interactions are 
present, at the working vortex density they do not much affect 
the dynamics leading to essentially the same finite size effects 
as for the EW line. At higher densities in-plane interactions 
should become the most relevant interaction of the problem. 
(As an example, a change in the scaling exponent was found 
with hard-core line interactions^) 



C. The vortex glass 

Here we analyze how the finite size effects reflect in the low 
temperature phase. We show data for T = 0.05 and y = 10~ 4 . 
We anticipate that aging effects appear at low temperature, but 
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FIG. 14: Dynamic wandering as a function of the plane index z. 
At this low temperature and for this disorder strength the system is 
aging. We show data for t w = 10 4 and each curve is labeled by At. A 
change in the slope of IV at short z is apparent. For this value of t w 
the longitudinal correlation length ^y(Af = 10 6 ) does not reach the 
system size L. 



FIG. 15: Normalized static structure factor Q(f,r w ) for At = t w = 
10 4 and different T and y given in the key. The value ko = 0.27 cor- 
responds to the first peak. The inset shows C%(t,t w ) for At = 10 4 and 
different waiting-times, t w = 10, 10 2 , 10 3 and 10 4 , with y= 10~ 4 
and T = 0.03. The structure factor does not vary significantly with 
waiting time while the system is aging. 



we postpone their discussion to Sec. lVIII In Fig.ll3lthe rough- 
ness and dynamic wandering evolution are shown. The char- 
acteristic aging t w -dependence is evident. The values of t \ and 
ti used in this figure are the same as in the EW case, Fig. |9] 
The first stationary single pancake regime is still present. Ag- 
ing develops in the single line regime but saturation seems 
to take place at the same value f^, suggesting that the satu- 
ration time t x does not strongly depend on temperature and 
disorder intensity for the studied parameter values. The dy- 
namic wandering also presents aging and the same saturation 
time is apparent. In Fig.^Jthe z-dependence of the dynamic 
wandering for t w = 10 4 is shown. Note that for this value of 
t w and the maximum At = 10 6 reached, the longitudinal cor- 
relation length does not reach the system size L, and then the 
Af = 10 6 curve shows a saturation regime for high z. Although 
the same trend as in the EW case is present, a difference is ob- 
served in the £j| | (Af ) < z regime. For instance, the curves with 
Af = 1 14305 and Af = 10 6 show a change in the slope. This is 
related to the fact that the roughness exponent £ of the Family- 
Vicsek scaling J45i changes with disordenSSS In the presence 
of disorder, at a given temperature there exists a characteristic 
length z c {T) separating two values of £ . For z < z c ■ ( T ) one has 
£ = £r, while for z > z c (T) one has £ = £d, with C,j < £ D , 
where £7 and Qj are the thermal (EW) and disorder scaling 
exponents, respectively. 

From the above considerations we choose to work in the 
time regime corresponding to both Af and t w being shorter 
than t x . This ensures that aging features do not overlap sat- 
uration effects and finite size effects should be irrelevant. For 
L = 50 this corresponds to times Af,f H . < 10 4 . 



VII. AGING IN THE VORTEX GLASS 

In this Section we describe aging features observed in the 
present simulations. As shown in Sec. IV Bl when entering 
the VG phase, for T < T g , is not possible to equilibrate the 
system, and correlation and response functions depend on the 
two times, f and f, v , showing aging. We start by describing 
the t w and fe-dependence of the correlation Q(f ,t w ). Then we 
discuss the aging characteristics of the pancake MSD B(t,t w ) 
and its associated integrated linear response %(t ,t w ). 



A. C k (t,t w ) 

Let us first present the ^-dependence of the normalized 
static structure factor for fixed Af and t w . In Fig. 1151 

we present data for Af = f„ = 10 4 and different T and y given 
in the key. At low temperature, in the VG phase, the curves 
present a well defined first maximum, followed by an observ- 
able second peak, but clearly not showing the presence of 
an ordered structure. At high temperature, well in the VL 
phase, the first peak is not as sharp as in the VG case, but it is 
clearly observed. For reference, the first peak corresponds to 
ko = 0.27. This is the value that we used in Fig. ^2 I n tne inset 
of Fig. El tne evolution of Q(f ,t w ) for fixed Af with waiting 
time is shown. Data correspond to y = 10~ 4 and T — 0.03. 
The structure factor is not strongly modified while the system 
is aging, i.e. while the relaxation time increases for increas- 
ing t w (see Fig. H6lbelow). This feature is also observed in 
simulations of structural glass formers using Lennard-Jones 
potentials i£ 

Figure IT6l shows the evolution of the correlation Q (f,f w ) 
at two temperatures, T — 0.2 > T g and T = 0.03 < T g , and for 
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FIG. 16: Correlation Q (f,? w ) for y = 10 and two temperatures, 
T — 0.2 and T = 0.03. Different waiting times are shown: circles, 
squares, diamonds and triangles correspond to t w = 10, 10 2 , 10 3 and 
10 4 , respectively. 



FIG. 18: Correlation C^(t,t w ) for fixed t w =10 and two temper- 
atures, T = 0.2 and T = 0.03 with y = 10~ 4 . Curves for different 
wave- vector values k are shown. 




10° io 1 io 2 io 3 io 4 io 5 



scaling. We give further support to this proposal below. 

The temperature dependence of Ck (t,t w ) for t w = 10 4 is 
shown in Fig.^] At the highest temperature the system is in 
equilibrium but for the remaining values it evolves out of equi- 
librium. The relaxation time decreases for decreasing temper- 
ature (at fixed t w ). 

FigurefTSlshows the ^-dependence of Q(f,f w ) forf w = 10 4 
and two temperatures, T = 0.2 and T = 0.03. At both temper- 
atures the system decorrelates faster with decreasing wave- 
vector k, meaning that the relaxation is much faster on small 
than on large length scales. This is exactly the same behav- 
ior observed in the Lennard-Jones glass-forming liquid out of 
equilibrium by Kob and Barrati^ 



FIG. 17: Correlation C ko (t,t w ) for fixed t w = IO 4 and y= 10~ 4 . 
Different curves correspond to different temperatures given in the 
figure. 



different waiting times. The disorder intensity is y= 10 . 
At high temperature the system decorrelates rapidly and the 
curves become stationary, i.e. they do not depend on t w . 
At low temperature a short time-difference quasi-equilibrium 
regime is present, but a clear f, v -dependence is developed at 
longer time-differences. The system decorrelates slower at 
longer waiting times, i.e. the system ages. For the longer t w 
value the curve does not clearly present a plateau as it usually 
does in other glassy systems. 36 An additive scaling scenario 
would require that the plateau develops at a longer t w . To 
exclude this possibility we should reach a longer t w while still 
being in the t w < t x regime. This demands to use larger system 
sizes L. We believe, however, that this is not necessary since 
the correlation actually decays according to a multiplicative 



B. B(t,t w ) 

In Fig. the time and temperature evolution of the pan- 
cake MSD 50 B(t,t w ) is shown for different waiting times and 
fixed disorder intensity y = 10~ 5 . B(t,t w ) is stationary for 
the high temperature value T = 0.2. For temperatures lower 
than T g w 0. 18 the aging behavior develops. As temperature is 
decreased the separation between curves with different wait- 
ing times increases, indicating that the system is deeper in 
the aging regime. Figure |20| shows the temperature depen- 
dence of B{t,t w ) for two waiting times, t w = 10 2 in Fig. l20r a) 
and t w — 10 4 in Fig. I20t b). After the initial single pancake 
regime, B ~ At, the sub-diffusive elastic line regime develops 
and B ~ At a . In the VL a = 1/2, while below some glass 
crossover temperature to be determined below, Cc(T) < 1/2. 
The scaling properties of B(t,t w ) will be analyzed in Section 
IVIIII The same behavior is observed using y = 10~ 4 . 



15 




t-t 

w 



FIG. 19: Pancake MSD B(t,t w ) at different temperatures. Circles, 
squares, diamonds and triangles correspond to t w = 10, 10 2 , 10 3 and 
10 4 , respectively. 



C X(t,t w ) 

In Fig.|^the dependence of the response function %{t,t w ) 
with t and t w is shown for y = 10~ 5 and T = 0.02 < T g . There 
is also an aging behavior, but not as severe as that of B(t,t w ) 
for the same parameters (Fig. \19\ lower panel). This fact is 
related to the type of scaling associated to the aging regime, 
which will be analyzed in Sec. lVIIII 



D. Comparison with aging of the single elastic line 

We present here results for the out-of-equilibrium dynamics 
of single elastic lines with disorder, i.e. without the in-plane 
vortex-vortex interaction. In Figs. 1221 and 1231 we present the 
evolution with time of B(t,t w ) and C(f,f w ) for the high and 
low temperature regimes, T — 0.2 and T = 0.05, respectively. 
At high temperature the dynamics is stationary, not showing 
any ^-dependence. At low temperature the system ages, as 
observed in the f w -dependence of both correlations in Fig. 1231 
These features were also observed in the related model of the 
directed polymer in random media, 71,90 and are qualitatively 
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FIG. 20: B(t,t w ) for different temperatures and for two waiting 
times: (a) t w = 10 2 and (b) t K = 10 4 . 




FIG. 21: Integrated response T%(t,t w ), where we include the tem- 
perature prefactor (kg = 1). The data correspond to 7= 10~ 5 and 
T = 0.02. 
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FIG. 22: Stationary dynamics in the high-T regime (T = 0.2) of 
the single line, i.e. without in-plane interactions, (a) B(t,t w ) and 
(b) Q (f,?vv) for y = 10~ 5 and different waiting times. The data for 
t w = 10 3 , 10 4 , 10 5 , and 10 6 collapse on a master curve. 



FIG. 23: Aging in the low-T regime (T = 0.05) of the single line 
case, i.e. without in-plane interactions, (a) B(t,t w ) and (b) Q (f,r w ) 
for y= 10~ 5 and different waiting times: t w = 10 3 , 10 4 , 10 5 and 10 6 
from left to right. 



the same as those presented in Sect. lVIll for the full interacting 
VG. This suggests that the out-of-equilibrium dynamics of the 
interacting vortex problem is dominated by the relaxation of 
the elastic line. We return to this important issue when ana- 
lyzing the FDT violation in the next Section. 



VIII. AGING SCALING AND FDT VIOLATION 
A. Multiplicative aging 

In this Section we study the two-times scaling of B(t,t w ) 
and #(f,f w ) in the aging regime as well as the violation of the 
FDT and the definition of the effective temperature. 

In Fig.|24]we show different scaling scenarii for B(t,t w ) at 
T = 0.02. On one hand, in Fig. 1241 a) an additive sub-aging 
scaling is put to the test by plotting B(t,t w ) vs. (f — t w )/tw, 
while the unsealed B(t,t w ) data is shown in the inset. On the 
other hand, in Fig. 1241 b) the multiplicative scaling of Eq. (|6) 
is used by plotting B = t" ■ B(t,t w ) vs. (t - t w )/t w . Although 
the sub-aging scaling with a low ji exponent seems appro- 



priate for long time-differences we prefer the multiplicative 
scaling for at least two important reasons. First, the bare data 
shown in the inset to panel (a) do not show a plateau at a con- 
stant B nor any clear trend to develop it at longer t w 's. This 
tends to disqualify the additive scaling for which a plateau 
is expected. Second, the multiplicative scaling allows for 
a better data collapse that includes the short time-difference 
sub-diffusive regime, only leaving apart the very short time- 
difference regime with normal diffusion. Thus, we find that 
our data is very satisfactorily described within a multiplicative 
aging scaling scenario. This scaling is similar to the one pro- 
posed for the low temperature behavior of the directed poly- 
mer in random media 71 and Sinai's diffusioni^^iZi Follow- 
ing Yoshinc 71 we tried the scaling form 

B(t,t w )=t%B{t) (46a) 



2k B Tx(t,t w )=Cx(t), (46b) 
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FIG. 24: (a) Test of a sub-aging additive scaling for the pancake 
MSD for y = 1(T 5 and T = 0.02. The bare data are included in 
the inset, (b) Multiplicative scaling for the same data as in (a). The 
values of the scaling exponents fl and a are quoted. A better data 
collapse is observed in the last case, including the short time regime. 
In the inset of (b), the scaled MSD and the integrated response func- 
tion are shown. Different waiting times collapse on the same master 
curve. The violation of the FDT at long scaled times t 2> 1 is clear. 



with t — 1 1 t w and the scaled variables B and % given by 

f cl (T)(?-i) a W f«i, 



B(t) 



C2 



(T)(f-l) a ^ 



f d(r)(?-i)«cr) 

X[t > \y(T)c 2 (T)(t-ir< 



t > 1 



f < 1 
t > 1 



(47a) 



(47b) 



where c\, and y are temperature-dependent coefficients (see 
the discussion in Sect. lVIIIBt . We assumed a linear depen- 
dence of x with B, but a non-trivial function % (B), as found in 
the scalar field or Edwards-Wilkinson equation in D = 2 and 
the spin-wave approximation to the 2D XY model, 64 - 65 cannot 
be discarded. 

In the inset to Fig.l24Tb) the scaled B and % at T = 0.02 are 
shown. The data for different t w fall on two master curves. B 
and x coincide for f<l, which means that FDT holds, while 
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FIG. 25: A single elastic line in a disordered environment: multi- 
plicative scaling of B(t,t w ) for the data shown in Fig[23ja). y = 10 -5 
and T = 0.05. The value of the scaling exponent a is given in the 
key. Data for different waiting times collapse, except for the station- 
ary normal diffusion regime at very short time-differences. 



for longer times f 3> 1, the value y(T) < 1 signals a violation 
of FDT, which we will analyze in Sec. lVIIICl 

In Fig.|25]the multiplicative scaling for the single line with- 
out in-plane interactions is shown. The curves correspond 
to y = 10~ 5 , T = 0.05 and different waiting times. Besides 
showing that this type of scaling is appropriate, this figure em- 
phasizes that the only regime not satisfying the scaling is the 
very short stationary diffusion regime, as suggested from the 
general picture in Fig. [2] 



B. Growing length 

Scaling arguments suggest that the averaged dynamics of 
an elastic line in a random environment should be determined 
by a single growing length scale, that separates equilibration 
at short length-scales from non-equilibrium at long length- 
scales. This is also, essentially, the picture that dictates dy- 
namic scaling in coarsening systems. 

The multiplicative scaling of Eg s. ( 146 at and d46bi implictly 
assume a growing length with a power-law dependence in 
time. However, the analytic determination of such a time- 
dependent length is very hard; it is known in a small number 
of problems such as ferromagnetic domain growth in clean 
systems with conserved and non-conserved order parameter." 
But even these presumably simple cases can be plagued with 
pre-asymptotic regimes, as shown for instance by Krzakala 
for the 2D Ising model with Kawasaki dynamicsilffi 

In the case of an elastic string in a random environment 
a phenomenological creep argument complemented with the 
assumption that typical barriers scale as L e leads to 



L(t) ~ L c 



u c U 



1/6 



(48) 
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asymptotically, i.e. in the limit of very long times and very 
large scales. L c is the Larkin length, U c its corresponding en- 
ergy scale and to a microscopic time scale. The exponent 9 
is usually further assumed to take the same value as for the 
free-energy fluctuations, 9 = 1/3. Alternatively, a similar ar- 
gument with logarithmically growing barriers imply a power 
law growth of the typical length scale L(t). 

Yoshino^Siffi and Kolton et studied the growing 

length-scale numerically in a lattice and a continuous model 
of a single line in a random environment, respectively. Both 
studies show that L(t) crosses over from a power-law to a 
logarithmic growth at a typical time-scale t* associated to 
a typical length scale L*. For instance, Kolton et al found 
t* ~ 10 4 and L* ~ 80 using lines with total lengths L = 256 
and L = 512. Similar values were found by Yoshino. 

In our study of vortices in interaction we were forced to 
use rather short total lengths, L = 50. Hence, we may suppose 
our dynamics falls in the first regime in which L(t) grows as 
a power law. Indeed, all our data can be satisfactorily scaled 
in a manner that is consistent with such a power law. We may 
expect, clearly, that the crossover in the growing length will 
induce a change in the scaling forms (|46aJ and J46bi possibly 
to a more general form 



B(t,t w ) ~ L^(t w )B a _ 
X{t,t w ) ~ lJ(t w )Xag 



( L(t) 
s \L(t w ) 
L(t) 



L{t w ) 



(49) 
(50) 



For times t < ?*, with a behavior L(t) ~ t a ^ and C = C one 
recovers Eqs. (146 a> and ( I46bl . For longer times a growing L(t) 
as in Eq. J48i is also possible. Testing the latest time-regime, 
however, goes beyond the scope of this study. The value of 
the exponent £ is not completely clear yet. Our numerical 
data support C = C but at longer times one might notice a 
deviation. 



C. Violation of FDT and effective temperature 

In order to test the violation of FDT, we use our results of 
Sec lVlllAI scaled with Eqs. ( 146 a> and ( I46bi . The parame- 
ter y(T) in Eqs. ( 146 at and d46bb measures the modification of 
the FDT, 2k B Tx(t,t w ) = y(T)B(t,t w ) (or % = y(T)B), and an 
effective temperature 37 can then defined by r e ff = T jy (see 
Eq.|42j. 

In Fig.[26]we show a parametric plot of the scaled variables, 
X(t) vs. B(t), for a temperature T < T g . If FDT is violated 
then the slope of the parametric plot should be different from 
unity. The parameters in Fig.[26]are y = 10~ 5 and T — 0.02, 
and different waiting times are shown. At very short rescaled 
times t all curves fall onto the FDT line, while for longer times 
the curves asymptotically change to a different slope y(T). 
In the inset, we show the short time regime for different t w 
values. Since the data with lower waiting time, t w = 10, is 
deeper in the out-of-equilibrium regime, we used the slope of 
X(B) with t w = 10 to compute the y(T) parameter signaling 
the FDT violation. In Fig.|^]we show x(B) with t w = 10 at 




FIG. 26: Parametric plot jtJ(B) in the VG (same data as in Fig. l24> . 
FDT holds at short rescaled-times while a violation of FDT with 
y < 1 appears at longer rescaled-times. The inset shows the short 
rescaled-time FDT regime for different t w values. 




FIG. 27: Parametric plot ^(fi) in the VG at different temperatures 
given in the figure. 



different temperatures; it is clear in the figure that the y(T) 
parameter decreases with decreasing temperature. These and 
similar curves are used to compute the y(T) values showed 
in Fig. I28f b). Exactly the same behavior is obtained for the 
single elastic line (not shown). 

In order to investigate the temperature dependence of dif- 
ferent quantities in the low temperature regime, we show the 
temperature variation of the scaling exponent a(T) and the 
y(T) parameter in Fig. [28] We use two values of the disor- 
der intensity, j— 10~ 5 and y = 10~ 4 . Also shown are data 
for 7 = 10~ 5 but for the case without in-plane interactions, 
quoted as single line. The values of the exponent and FDT 
violation factor at high temperature are a = 1/2, correspond- 
ing to the single elastic line without disorder limit (or VL 
phase), and y = 1, corresponding to the equilibrium FDT. In 
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Fig. l28f a) we see that a(T) depends weakly on T and that it 
decreases with increasing disorder strength within the glassy 
regime. It is also observed that the in-plane interactions tend 
to decrease the value of a. The parameter y(T) measuring 
the violation of FDT rapidly increases with increasing tem- 
perature. In Fig. I28t b) we show that y(T) = T /T e g is well 
described by a linear form, implying an effective temperature 
T e g that is temperature-independent. However, a more com- 
plicated temperature-dependent effective temperature cannot 
be ruled out. We do not show the linear fit corresponding to 
y = 10~ 5 without in-plane interactions. From the linear fit 
of y(T) one obtains the effective temperature values quoted 
in the figure, r eff = 0.175 and T eS = 0.203 for y = 10~ 5 and 
y = 10~ 4 , respectively. 

Another feature of the scaling form ( 147 \ is the temperature 
dependence of the c\(T) and C2(T) coefficients. These coef- 
ficients measure the separation between the quasi-equilibrium 
regime and the aging regime, and hence they should converge 
to the same value in the equilibrium high-temperature regime. 
In Fig.|29]these functions are shown for the same parameters 
used in Fig. [28] The T e ff values are also given in the figure. 
For each set of parameters c\ (T) and C2(T) are closer to each 
other when approaching T e f[ and the high-temperature regime. 
This is another confirmation of the scaling form j47l >. 

Noticeably, the obtained value of T e ff = 0.175 is very close 
to the crossover temperature T g w 0.18 below which a dy- 
namic arrest is observed and the system can not be equili- 
brated, as found in Sec IV Bl A similar result is observed in 



structural glasses: T e g ss T g (as in a random energy model sce- 
nario —). 



IX. DISCUSSION 

In this paper we studied numerically the dynamics of a 
high-temperature superconductor model. We briefly analyzed 
the stationary dynamics in the VL and we focused on the out- 
of-equilibrium dynamics of the VG. Our aim was to learn 
about the nature of the VG and its dynamical properties from 
the relaxation of different two-times correlation functions. 
During this study a systematic comparison with the dynam- 
ics of other glassy systems was performed. 

A key feature in problems involving elastic lines, as the 
one we treated here, that should be highlighted is the rele- 
vance of finite size effects. By comparing the evolution of the 
roughness and dynamic wandering in the vortex system to the 
ones in the EW equation for growing interfaces 89 we identi- 
fied three dynamic regimes present at all interesting tempera- 
tures: (z) a very short time-difference normal diffusion regime 
without effective elastic interactions; (ii) a sub-diffusion in- 
termediate time-difference regime characterized by a grow- 
ing longitudinal correlation length; (Hi) a long time-difference 
regime where the correlation length has reached the system 
size and normal diffusion simply reflects the center of mass 
diffusion. In order to separate aging effects from finite size 
effects we constrained the remaining study to regimes (z) and 
(zz), i.e. before roughness saturation. 

Strikingly, the behavior of the density-density correlation 
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FIG. 28: (a) The scaling exponent u(T) and (b) the y(T) parameter 
measuring the FDT violation in the VG with two disorder intensities, 
y= 1(T 5 and y= 1(T 4 . Data for a single elastic line in a random 
environment with y— 10~ 5 are also shown. The effective tempera- 
ture r e ff obtained from linear fits y(T) = T/T e g is given for the cases 
y= 1(T 5 and 7= 1(T 4 . 




FIG. 29: Evolution of the coefficients Q (T) and C2(T) with temper- 
ature for the same parameters as in Fig. 1281 They tend to the same 
value with increasing temperature, as suggested by the scaling form 
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function Q(f,f H ) resembles in many aspects the one observed 
in simulations of Lennard- Jones glass formers. 62 In particu- 
lar, the wave-vector dependence of the correlation follows the 
same general trend. The main difference is that the density- 
density correlation does not develop a well defined plateau for 
long waiting times in the VG, while it seems to do in glass- 
forming liquids. This is not due to the use of extremely short 
waiting times as could be expected in an additive aging scal- 
ing, but it is in the core of the multiplicative aging scaling we 
found in the present simulations. 

We showed that the two-times evolution of the pancake 
mean-squared-displacementB(f ,t w ) is very well described by 
the multiplicative scaling diffusive-aging scenario earlier pro- 
posed in the numerical study of the out-of-equilibrium dynam- 
ics of the directed polymer in random media model, 71 found 
analytically in the massless scalar field in D = 1 (related to the 
Edwards-Wilkinson surface), and obtained numerically and 
analytically in Sinai diffusion. 75 ' 76 ' 77 It is interesting to stress 
that the model of a D-dimensional manifold embedded in an 
Af-dimensional space under the effect of a combined disor- 
dered and harmonically confining potential - a toy model that 
is usually used to model vortex systems - fails to capture the 
multiplicative scaling in the large limit. 73 - 74 This model has, 
instead, a well defined transition at a finite temperature below 
which correlations and responses scale in an additive way with 
a stable plateau. 

We studied the violations of the FDT out of equilib- 
rium. We constructed the parametric plot of linear inte- 
grated response against displacement after having eliminated 
the power of the waiting-time, t", that appears multiplying 
these quantities. The resulting plot is linear, at least within the 
accuracy of out data, and it allowed us to identify an effective 
temperature as the slope of the linear plot. The same behavior 
was found in the lattice model of the directed polymer in a 
random environment. 69 - 70 The relation between integrated re- 
sponse and displacement, once the factors t " have been taken 
into account, might involve though a non-trivial function of 
the displacement itself, not visible within the accuracy of our 
data, as has been found, for instance, in the 2D scalar field 
and XY model. 64 - 65 The FDT violation in the ' mean-field' 
model of a relaxing D-dimensional manifold in an infinite- 
dimensional embedding space 73 - 74 with a short-ranged disor- 
der potential has this type of violation of FDT though it is not 
necessary in this case to eliminate additional t„ factors. 

We performed the analysis of displacement, correlations 
and responses at several low temperatures and we analysed 
the temperature dependence of the exponent 0f(T) as well as 
the coefficients c\(T) and C2(T) that characterize the scaling 
function. We also followed the evolution with temperature of 
the parameter signaling the FDT violation, y(T). We consid- 
ered two intensities of the pinning disorder and we compared 
the data to the a priori simpler case of a single line without 
in-plane interactions. All the parameters, y(T), 0C(T), c\(T) 
and C2(T) show a clear trend to reach the high-temperature 
behavior: y(T) = 1, a(T) = 1/2, and c x (T) = c 2 (T) when 
approaching the VL from the glassy regime. From these out- 
of-equilibrium measurements the crossover temperature was 
found to be very similar to the value of the effective tempera- 



ture, r e ff, obtained from the fit of the FDT violation at very 
long time-differences to a linear form with y(T) = T/T e ff. 
This result is similar to what happens in the random-energy 
model in which the effective temperature takes the transition 
value in the full low-temperature phase. 

To further test the meaning of the crossover temperature and 
the effective temperature T e g, we studied the dynamic arrest 
approaching the VG from the VL. By fitting a mean-squared- 
displacement correlation, taking into account sub-linear diffu- 
sion, we determined a relaxation time t r . The so-obtained re- 
laxation time increases with temperature as t r ~ T~ 2 , at rather 
high temperatures, as expected in disordered free dynamics. 
However, below a characteristic temperature a sudden further 
increase of the relaxation time f, was found. We identified the 
crossover temperature between normal relaxation and rapidly 
growing relaxation with a glass temperature T g . When com- 
paring with the effective temperature we found that both tem- 
peratures are of the same order, a fact already observed in 
other glassy systems. This suggest a freezing of the slow de- 
grees of freedom of the system around the crossover tempera- 
ture. 

One outstanding result of our simulations is that all the out- 
of-equilibrium properties are mainly dictated by the elastic 
line energy, i.e. they are all observable in the single line with- 
out in-plane interactions. This means that for the parameters 
used in our simulations the two relevant competing energies 
are the elastic line and pinning energies. This fact points out 
that the relevant dynamic properties are also observed in sim- 
plified models of elastic lines in disordered media as directed 
polymers2ii22i2P. or continuous models similar to that studied 
here42ii21 It is worth mentioning that this cannot be the case 
at very high magnetic fields, i.e. very high vortex density, 
since at some point the in-plane interaction should become the 
relevant energy scale. In this condition strong excluded vol- 
ume effects, such as the "cage effect" in supercooled liquids, 
should be appreciable. A simplified model in this direction 
was considered by Petaja et al., who considered elastic lines 
with hard-core interactions in random environments. 97 

From our findings several extensions could be envisaged. 
Concerning the possible experimental observation of the out- 
of-equilibrium dynamics of the VG phase, new experiments 
monitoring magnetic relaxation 44 - 45 would be most welcome. 
To measure aging, one possibility is to perform transport ex- 
periments following a protocol similar to the one used by 
Ovadyahu and coworkers to study the electron glass regime 
in Anderson insulators .iffiiiS^ Tests of FDT would require si- 
multaneous measurements of noise (magnetic and/or voltage 
noise) and time dependent response (relaxation of magnetiza- 
tion and/or resistivity); similar tests have been performed in 
other systems piSiiSSJS 6 . Another possibility would be to test 
the aging dynamics and the FDT violation in vortex shear ex- 
periments. 107 - 108 It is well known that in shear-thinning sys- 
tems the time scale introduced by the shear rate basically 
replaces the waiting time; at lower shear rates the relax- 
ation is longer, i.e. the system "ages" with decreasing shear 
rate> i09 i iio i ui 4 ii2 Qj t j]j an0 ]- ner possibility is to measure effec- 
tive temperatures using a tracer embedded in the vortex sys- 
tems. Finally, the out-of-equilibrium dynamics of the center 
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of mass diffusion could be experimentally tested. 

It would also be interesting to study the out-of-equilibrium 
relaxation of the BG and compare the outcome to the results 
obtained in this work for the VG. 
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